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The standard approach to initial data for both analytic and numerical computations of black hole 
collisions has been to use conformally-flat initial geometry. Among other advantages, this choice 
allows the simple superposition of holes with arbitrary mass, location and spin. The conformally flat 
restriction, however, is inappropriate to the study of Kerr holes, for which the standard constant-time 
slice is not conformally flat. Here we point out that for axisymmetric arrangements of rotating holes, 
a nonconformally flat form of the 3-geometry can be chosen which aUows fairly simple superposition 
of Kerr holes with arbitrary mass and spin. We present initial data solutions representing locally 
Kerr holes at large separation, and representing rotating holes close enough so that outside a common 
horizon the spacetime geometry is a perturbation of a single Kerr hole. 
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I. INTRODUCTION AND OVERVIEW 



The study of black hole collisions has been of recent interest in connection with the development of gravitational 
wave detectors [Q], and as a proving ground for numerical relativity In addition, it will give unprecedented insights 
into the nonlinear workings of relativistic gravitation in strong field situations. 

The approach in numerical work has usually been to use the prescription of Bowcn and York to generate initial 
data representing the holes. In this prescription, the initial value problem is greatly simplified by the assumption 
that the initial geometry is conformally flat. This approach, though, has some disadvantages for studying collisions 
with net angular momentum, since the final hole formed will be a Kerr hole, and — at least for standard slicings 
— the spatial geometry of a Kerr hole is not conformally flat. This means that computations of radiation from 
coalescence will contain early radiation due to the relaxation of Bowen-York holes to Kerr holes, as well as radiation 
more directly attributable to the coalescence itself. Still, this represents an astrophysical imperfection rather than a 
major obstacle to numerical work. For analytic work, however, it is useful to consider initial configurations close to 
the final stationary hole so that the initial configuration, and the evolution, can be treated with perturbation theory 
This technique is inherent in the methods used for extraction of radiation ||^ in full numerical relativity. First 
and second-order perturbation theory have also been remarkably successful in the "close limit approximation," the 
analysis of collisions starting from small separations [|6|-^ . 

In this paper we present an approach to the generation of initial data that fits the needs of perturbative analysis. 
Since this approach has essentially the same computational simplicity as that in the conformally flat approach, it is 
our hope that these solutions will be considered as starting points for numerical relativity evolution. In recent work, 
comparisons @Jl^ between analytic and numerical results have been very helpful to both efforts. 

Within the context of a more general formalism, we will present two families of initial data solutions. The first 
family of solutions will represent an axisymmetric superposition of two equal mass Kerr holes, parameterized by the 
mass of each hole and the spin of each hole, and a measure of the separation between them. The solution will have 
the crucial feature that the holes have local "Kerricity." That is, in a region around each hole small compared to the 
spacing between the holes, the 3-geometry and extrinsic curvature will have the form of the Kerr initial data, aside 
from corrections that vanish in the limit of infinite separation. The specific solution given serves to illustrate a method 
of superposition which can be generalized in an obvious way to a large number of axisymmetric initial configurations. 
In particular, it can be used to construct initial data for unequal (different mass, different spin) Kerr holes, for more 
than two holes, for a collection of Kerr and Bowen-York spinning holes, or for a collection of Bowen-York holes 
(although in this last case the method involves only conformally flat geometries and is equivalent to the Bowen-York 
method) . For any solution constructed by this method the hole will locally have the same form as if it were isolated. 

A very important feature of these solutions is their close limit. If we reduce the separation between two holes, an 
initial horizon will surround both holes, but in the limit of small separation the solution outside that horizon is not 
a single Kerr hole. This is very different from what happens for nonrotating holes, and this difference is one of the 
major points we wish to make in this paper. The fact that the Kerr geometry is not the "natural" close limit for a 
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sequence of initial data sets for two locally Kerr holes is not a flaw in our superposition method, but is to be expected 
for general initial data sets. 

The fact that the close limit of rotating holes is not a Kerr hole is surprising because it is so very different from 
the Schwarzschild case. For initial value solutions for two nonspinning holes, such as the solutions of Misner |l6| or of 
Brill-Lindquist [ p^f , if the separation between the holes is reduced one expects an all-encompassing horizon to exist on 
the initial hypersurface, and for the initial value solution outside that initial horizon to approach the Schwarzschild 
geometry. And this is just what happens [^|8j. In the mathematics, as in intuition, the only source of nonsphericity 
is the separation of the initial holes, and as that separation is reduced, the solution outside the event horizon becomes 
spherical, and hence becomes Schwarzschild initial data. 

Crucial to the simplicity of this process is the fact that the Schwarzschild initial data is singled out by the requirement 
of sphericity. The Kerr initial solution, by contrast, requires much more: a set of relationships involving the initial 
extrinsic curvature as well as the initial 3-geometry. If these relationships are satisfied for the individual holes, there 
is no reason for them to be satisfied by the superposition of two holes. One can see a rough quantitative sign of this 
in the following. For a rotating hole of mass M and angular momentum J, the quadrupole moment in the 3-geometry 
has a magnitude of order J'^/M. If we add two such holes we expect to infer from the distant geometry that we have 
approximately twice the mass, twice the the angular momentum, and (from the part of goo) twice the quadrupole 
moment. The quadruple moment will therefore be too small by a factor of two for the result to be a Kerr hole. In 
a manner of speaking, then, what makes two Kerr holes different from a single hole is not simply the separation, it 
is the more basic fact that there are two of them. There is a completely different point of view from which this can 
be understood. The Kerr solution is singled out as a spacetime geometry with a stationary horizon. The initial value 
solution — the 3-geometry and extrinsic curvature — are not special as initial value solutions. By contrast, the initial 
value solutions for Schwarzschild (spherical 3-geometry, zero extrinsic curvature) is special. 

We are pointing out here that the Kerr close limit should not be expected from the superposition of two Kerr holes. 
This should not be confused with the fact that the Kerr geometry is the natural late time result for the evolution 
of two holes. Due to the evolution, the geometry and extrinsic curvature on a late hypersurface, outside the final 
horizon, will have approximately Kerr form. For that reason, it is useful to have solutions available with a Kerr close 
limit. 

In Sec. Ill below we construct such a solution. This solution is somewhat contrived, but any other multihole 
solution with a Kerr close limit will also be contrived (except of course, the solution produced by evolution). Our 
particular contrived solution has the unappealing feature that there is a singularity where there should not be a 
physical reason for a singularity. Since this contrived solution is meant to be used only for small separations, in which 
case this singularity will be well inside the horizon surrounding the two initial holes. In judging the importance of this 
singularity, it is necessary to keep in mind the use to which the close-limit solution is to be applied. For this initial 
value solution (or any other initial value solution with a Kerr close limit) it is possible to use close-limit perturbation 
theory to find the outgoing radiation generated in the evolution of the initial data. If these initial data are used for 
numerical relativity, very useful comparisons with analytic results are possible. For such comparisons the unwanted 
singularity is of concern only if it causes problems in the numerical method. For numerical relativity computations 
that use horizon excision the singularity should be of no concern whatever. 

An alternative approach has recently been given by Baker and Puzio [p^ to the construction of initial data sets for 
axially symmetric non-conformally flat, initial data. That method starts with a specification of the 3-geometry and 
solves for partial information about the extrinsic curvature. The method has the advantage, in principle, of allowing 
generalization to a broader class of axisymmetric solutions, in particular to rotating holes which are initially moving. 
The method has the disadvantage that it is more difflcult to implement, and has not yet been tested numerically. 

The remainder of the paper is organized as follows: In Sec. II we present the basic method of superposition, and 
the numerical solution for two identical "locally Kerr" holes, including a numerical demonstration of local Kerricity. 
In Sec. Ill we show that the locally Kerr solutions of Sec. II do not have a Kerr close limit, and we introduce solutions 
that do have a Kerr close limit. A discussion and conclusions are presented in Sec. IV. 



II. MATHEMATICAL APPROACH 

We start with the Brill- wave |^ form of the metric, as studied by Brandt and Seidel |2l| - [2^ . In place of their 
radial coordinate rj which is appropriate to the description of both sides of a throat, we use the alternative coordinate 
f, related to r/ by drj = df/f. In terms of this coordinate, and the Brandt-Seidel notation, the form of the metric is 

ds^ = (df^ + f^de^) + f2 sin^ d dcj)"^] , (1) 
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where the q and the conformal function $ are functions of r and 9. The extrinsic curvature has only rcj) and 6(j) 
components, and in the Brill-Brandt-Seidel notation those are denoted in terms of He and Hp, conformally related 
functions of f , 9 as 



'<^>-^ He sin 



Ked, = f-^^^"^ Hp Hint 



(2) 



The initial value equations of Einstein's theory impose no constraints on q, but the momentum constraint requires 
that the Hs satisfy 



^df{HE) sin^ 9 + deiHp sin^ 



= . 



(3) 



The fact that in this form the momentum constraint is linear, and does not depend on the spatial geometry (i.e., is 
independent of both q and is what ultimately will allow the fairly simple superposition to be described below. 
Once q, and Hs solving (0), are chosen, one finds the conformal function $ by solving the Hamiltonian constraint 



/:ham($) = \hI sin^ 9 + Hi]/ (4f6) , 



(4) 



where 



9^ 

sin 9 89 V''''' 89 



4f2 



dr \ dr'^ 



(5) 



Note that the terms in curly brackets have the form of the flat space Laplacian. 

For a Kerr hole of mass M and angular momentum parameter (angular momentum per mass) a, the forms of the 
functions arc 



cos^ 9 



„2 I 2Ma^rsin'^ 
~r — T"; — 7~ 



r^+a-^ cos^ 



(6) 



Hex 



aM [(r 



2)(r2 + a2 cos2 6») +2r2(r2 + 0^)] 



(7) 



H 



-2a^Mr { - 2Mr + a^) cos 9 sin^ I 



FK 



(8) 



K 



2Ma^r sin^ 9 



(9) 



Here we have used the index "if" to denote the Kerr form, and r is a rescaling of r given by 



r(f, M,a) = f [ I + 



M 



2r 



1 



M -a 
2f 



(10) 



The solution of (H) depends not only on how we choose q, He and Hp, but also on the boundary conditions for <i> 
that are specified. To recover the Kerr solution, we must not only specify that q, Hp and Hp have the forms given 
in (|6|)-(^), we must also specify boundary conditions that $ ^ 1 at r — > oo, and that $ has a singularity at f = of 
the form $ 



VA/2-a2/(2f). 

To specify more general boundary conditions for (|4|), it is useful to write 



(11) 



in which ^rcg, is regular in the finite f, 9, (p space, and $rcg — > at f — > cx3. The physical information specifying the 
nature of the solution to (||) that we are seeking is put into a specific choice for $sing- Equation (^ then takes the 
form 
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-'ham 



(*sing) - («'sing + *ro 



/(4r-6) . 



(12) 



This equation can be solved iteratively for $rcg, and can be viewed as a linear partial differential equation, with a 
known source, at each step of the iteration. The linear difference equations arising at each iterative step were solved 
with an explicit solver from the LAPACK package, and as a check we also solved with a successive overrelaxation 
routine. A compactified radial coordinate r,. = f/{f + c), reaching a finite value at f ^ cxi, was used, and it was found 
that values on the order of c = 1 tended to give the best accuracy. 

It is instructive to "find" the Kerr geometry with the method just outlined. We start by taking q, He and Hp to 
be the Kerr functions given in (j6|) — Next, a choice must be made for <i>sing- A choice of the form 



$sing = 1 



2f 



(13) 



will lead to a single throat solution. If, in addition, we specify that n — — a^, we are choosing 3>sing to have 

the singular behavior of ^k, and the solution for $ should be the Kerr conformal factor. The solution for any other 
nonzero choice of k may be considered to be a distorted Kerr hole. [If n is taken to vanish, the right hand side of ( |T^ ) 
diverges, and the solution for $ turns out to give a nonstandard topology, of a type discussed in the next section.] As 
a test of our numerical code, we have solved ( p^ ) for Kerr inputs qx, Hek, Hpx with the singular solution in ( p^ ) and 



the Kerr choice 



Numerical solutions were computed for the case M = 1, a = 0.9, and were found to 



be in excellent agreement with the analytic form in (|9|). For a grid with 200 radial and 80 angular divisions, the 
fractional error was of order 3 x 10~^. It should be noted that ^ 1 + M/(2f) at large f, whereas we have made 
the choice ^sing ^ 1 + V — a? / (2f). The difference in the asymptotic behaviors is supplied by $rog which must 
fall off as $reg {M — \/ — a^)/(2f). We have numerically verified, with an accuracy of 0.1%, that the coefficient 



of the 1/f term in the computation of ^k, is M/2. 

As the first step in superposing Kerr holes, we take a Kerr hole, as described above, for mass Mi and spin parameter 
Oi, and in the following way we shift the hole to a coordinate position f — zi,9 — 0: We introduce coordinates fi, ^i, 4>, 
related to our fundamental coordinates f , 9, (j), by 



ri cos Oi — r cos ( 



Ti sm Oi = r sm f 



(14) 



In the new coordinates the metric has the same form as in (|l|) so a solution of the initial value equations is given by 
® - (|l^) , with f , 6 replaced by fi , , and with the mass and angular momentum parameters taken to be Mi , ai . 
We now reexpress this shifted solution in terms of the fundamental coordinates f , 9; we denote the resulting solution, 
with a subscript 1, as ^i, qi, Hpi, Hpi. The first two functions transform rather simply, and are given by 



qi{r,9) = qK{ri,9i;Mi,ai) 



$1 = $x(ri,6'i;Mi,oi) , 



(15) 



where q^, are understood to be the functional forms of coordinates and parameters given in (|^) and (|^). In other 
words, these functions are given explicitly by 



-2qi 



rf + COS^ 9i 



(16) 



and 



2Mia?rsin^ 9i 



'•f + af cos^ 9i 



In these expressions, the symbol ri is to be interpreted as 

Ah + ai 



ri 



ri 1 + 



2fi 



1 



Ml - a I 
2fi 



(17) 



(18) 



and ri, 9i are taken to be the functions of r, 6* given by (|l4|). 

To find the shifted extrinsic curvature we simply transform the components of the tensor Kij and reinterpret the 
results in terms of the conformally related curvature quantities Hpi, Hpi. This gives us 



'■^i^Hei sin^ 9 = -^f^H],^HEK sin^ 9i + -^f^^^-^^HpK sm9i 



(19) 



4 



When these are combined with the second relation in (|l5|) we find 

2 



Hei sin^ 61 = — 



dri ^ 2 - ^ 

a?' or 



(20) 



(21) 



sin 6 



ri 



HpK sin 01 



(22) 



In these equations, the expressions Hek and HpK are understood to be the functional forms for the Kerr geometry. 
That is, they are given by (0) and (^) with f,0,M,a replaced by fi, 6*1, Mi, oi. We have used Maple to verify that 
the solutions in ([2l|),(p^ satisfy the momentum constraint (^). 

As a check on our numerical program we have used the shifted solution of (16), (^, and ^ 
the singular solution 



in (12) along with 



$sing = 1 



\Jr'^ — 2z\r cos Q + z\ 



(23) 



Equation ( |l2|) was solved numerically for <i>iog- The result for <i>iog is shown in Fig. 0, and is compared with the 
regular part of $i, the analytic solution for a single shifted hole given in (|7|), for a hole with M — 0.5, J — 0.225, 
and located at zq = 0.1. Part (a) shows the numerically computed result for <i>""g™ as a function of r and d\ an inset 
shows some detail of "fine structure" in the solution. The fractional difference between the computed and the analytic 
solutions is plotted in part (b). These computations used a numerical grid with 1500 radial and 800 angular grid 
divisions. The results, both numerical and analytic, show that *i>rcg is maximal near the location of the throat, and 
— as shown in the inset — has some fine scale structure where the symmetry axis passes through the throat. The 
fractional difference between the computed and the analytic <&reg was found to exhibit second order convergence and, 
at the finest grid spacing, was less than 0.5% even near the axis, where $reg changes on a small length scale. 

In a manner similar to that just described, we can find the functions ^2,^2, He2, Hp2, describing a hole of mass 
M2 and spin parameter 02, shifted to Z2- We need only replace all "1" subscripts in our first shifted solution by "2." 
Since the initial value equations are not linear we cannot linearly superpose the two shifted solutions, but we can 
superpose the q, He and Hp functions. That is, we can take 



q = qi + q2 



He = H 



El 



H 



E2 



Hp — iJ, 



Fl 



H 



F2 



(24) 



The first of these superpositions is immediate, since g is a "free choice" that is not required to satisfy any constraint. 
The H functions must satisfy the momentum constraint (|^), but that constraint is linear and independent of the 
3-geometry. Since the first and the second shifted solution individually satisfy the linear, geometry-independent 
momentum constraint, their superposition will satisfy it. 

It remains to solve for the conformal factor <I>, and this requires specifying the singularity structure. We do this by 
taking 



- 1 = ($singl - 1) + ($sing2 - 1) 



(25) 



where $singij *I'sing2 are the singular parts of $1, <1'2- This is equivalent to taking a Newtonian- like solution 



sing 



1 

1 + ::: 



2 \Jr'^ ~ 2z\r cos B + z\ 2 — 2z2r cos 9 + z\ 



(26) 



The last step is to solve ( |12|) numerically for <i>rcg- 

With the numerical solution of the Hamiltonian constraint ( p^ ) , we have an initial value solution for two holes, 
parameterized by the masses, spins, and locations of the holes. A crucial feature of this solution is that the spacetime 
is "locally Kerr" near the holes. Let us consider hole 1, and let us suppose that z\ — Z2 is large compared to M\. 
Kerr-like coordinates can be introduced in terms of which the 3-geometry and extrinsic geometry components have 
approximately the form of a Kerr hole with parameters A^i,ai. At a radial coordinate distance r from hole 1, the 
fractional deviation from Kerr form will be of order r^/jzi — Z2p, except for the influence of hole 2 in the conformal 
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factor. This correction will be of order r/|zi — 02 1- The proof of these statements is straightforward, but tedious, 
and the details are not given here. Instead, we present a numerical demonstration of local Kerricity for the case of 
identical holes (Mi = M2 = 0.5, Ji = J2 — 0.225) symmetrically placed about the coordinate origin {zi = —Z2 = zq). 

In Fig. ^, part (a) shows the computed regular part of <&, as a function of f and 9. Part (b) shows 

A = ($"™' — $reg)/^'reg, the relative difference between the computational result and the analytic function $rog, of 
( p7| ) for a single isolated hole of M = 0.5, J = 0.225. Since the (two-throat) singular part for the two hole solution 
differs from the (one-throat) singular part for an isolated solution, the difference shown in (b) is not the same as 
the fractional difference between the and 'I'rog, the regular part of the shifted Kerr solution. Rather, part (b) 

indicates the size of ($;™' - *rcg) = *"cg" - *rcg + \/M^ - a\l(2f2) with fa = - 2z2f cos 6* -h z|. In Fig. | 
the separation parameter is zq = 2.0, and the computed hole exhibits local "Kerricity," i.e., the deviation from the 
Kerr geometry is small, less than 4% in the neighborhood of the throat. Figure || shows the analogous results for a 
separation parameter zq = 0.25. Here, the deviation from the Kerr geometry is significantly larger, on the order of 
12%. Yet smaller values of zq produce yet larger deviations from the Kerr geometry. 

For definiteness we have described the axisymmetric superposition of two equal mass, equal spin, Kerr holes, but 
it should be clear that the method of superposition is more widely applicable to the axisymmetric superposition of 
solutions. The general principle is to put the initial value solutions to be added into the form (|l]). The individual 
metric q functions and extrinsic curvature H functions are then superposed as in (p^, and the momentum constraint 
(^ is guaranteed to be satisfied. The last step is to choose the singularities, and (numerically) to solve the Hamiltonian 
constraint for $rcg- This can be done, for an arbitrary axisymmetric array of Kerr holes, but the method is not limited 
to Kerr holes. Any solution that can be put in the form (|l|) can be superposed. In particular, the rotating holes of the 
Bowen-York 1^ formalism can be used as building blocks. For these conformally flat solutions, 9 = 0, and superposing 
only Bowen-York holes will lead to a conformally flat initial value solution; in this case our method is simply a special 
case of the more general superposition that is valid in the Bowen-York formalism. More interesting is the possibility 
of superposing arbitrary mixtures of Bowen-York and Kerr holes, or Kerr holes with distorted Kerr holes. 



III. CLOSE LIMIT SOLUTIONS 



For solutions constructed by superposition as described above, the close limit is not a single Kerr hole. In the two 
Kerr hole solution, for example, if we let \zi — Z2I 0, the initial solution will have an all encompassing horizon 
surrounding both the throat at zi and at Z2- Outside this horizon the initial value solution will not approach the Kerr 
solution. We illustrate this in Fig. ^. In part (a) we plot the computed conformal factor $ for the superposition of two 
identical holes, with masses Mi — M2 — 0.5 and angular momenta Ji = J2 — 0.225, placed at zi — —Z2 — zq = 0.02. 
In part (b) we plot A, the relative difference between the computed solution and a Kerr solution. The comparison 
Kerr geometry is chosen to have angular momentum J — Ji + J2 = 0.45. For any other choice the extrinsic Kerr 
geometry would not agree with that for the two hole solution at large f. By looking at the large f behavior of the 
computed $ it was determined that the two hole solution had a mass of Af = 0.824. (This is less than Mi + A/2 = 1 
as should be expected.) If we were to compare our solution for <f> with a Kerr hole of any other mass, there would be 
unacceptable differences at large f . 

In Fig. ^, the left (small radius) boundary of the surface is at the approximate location of the initial horizon. As 
one would expect, the plot in (b) shows that the maximum deviation from the Kerr geometry (around 8%) occurs on 
the horizon. The crucial feature of these plots is that they do not change significantly as zq — > 0. The close limit of 
the two holes is not a single hole; the deviation from the Kerr geometry exhibited in part (b) does not decrease with 
decreasing separation. This is in clear contrast to the nonrotating case, in which the deviation outside the horizon, 
decreases monotonically with the separation of the holes. 

There are two very different reasons for wanting a solution with a single-Kerr close limit. The first is that the 
evolution of a binary system will result in a single Kerr hole as the evolved close limit, so we might want an initial 
value solution with this feature. (This will be discussed further in the next section.) Of more practical interest, if we 
are to do close-limit perturbation theory, we must have a simple stationary solution as the close limit. 

It is, in fact, very simple to construct a superposition with the desired close limit. We can take the metric function 
q and the conformally related extrinsic curvature functions H to be precisely those of a single Kerr hole of mass and 
spin M, a, i.e., the functions given by (^) — (^) and (|o|). We could then try to impose the topology of two throats 
by taking $sing to have the two-throat form of (p6|). But in this case, the righthand side of ([l^ ) would have a 
singularity at f ^ 0, if ^sing + $rog were finite there, since Hek does not vanish at f ^ 0. In fact, the "nonsingular" 
function <i>i.eg we compute diverges as 1/Vf as f ^ 0. With the resulting $^ ^ behavior of the conformal factor, 
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the topology of the metric is Ti} x near r ^ 0. It is possible in principle for us to make $rog truly regular by adding 
a term of the form h{0)/y/f to $sing- The angular function h{0) is the solution of the nonlinear ordinary differential 
equation that results when the Kerr forms of He, Hp and q, and the ansatz $ = h[9)/\/f are put into (^), and the 
f — > limit is taken. Since this is numerically inconvenient, we have chosen to keep the 1/Vr singularity as a part of 
$reg, and to use the boundary condition for ( [l^ ) that \/f<i>rog is regular at f ^ 0. 

An alternative is for us to choose to add a term of the form n/f to <I>sing- For any value of k this will cause the 
extrinsic curvature term on the righthand side of (^2|) to vanish at f — + 0. Of course, by adding such a term to $singj 
we are adding another asymptotically flat region, i.e., another throat. Our solution, then can be taken to represent 
an "extra" black hole. We could minimize the intrusiveness of this unwanted additional throat by choosin g k , the 
parameter that governs its size, to be arbitrary small. In this case, at small f the extrinsic curvature term in ( |l2[) will 
become very large, of order J^/k^ around f ^ k. 

If we are to use the above choices of topology, either the $ ^ l/v^, or <I> 1/f, at f — > 0, we no longer have the 
two throat topology that represents two holes. But if this superposition is to be made only in the case the two holes 
are quite close, the additional central topological feature will be deep inside an all-encompassing horizon, and will be 
irrelevant to close-limit analysis or for numerical relativity if horizon excision is used. As we shall argue further in the 
next section, this means that the unwanted central topological feature in no way makes our close-limit superposition 
less useful. In this section we give numerical results for close-limit superpositions, and further details of how they can 
be used for numerical and analytic work. 

In our close limit superposition to form a single Kerr hole of mass M and spin a, in addition to choosing the Kerr 



forms for q^Hp and Hp, we must also choose $sing to have the form of (26), with the numerators in the second and 
third terms replaced by 1/2 ^/M"^ — a^. It is only for this choice of singular structure that the mass contained in $sing, 
combined with the mass coming from the solution of (p^, will be the mass M in the zq — > limit. 

In Fig. I we illustrate the solution to ( p^ ) for the problem just described. For two close holes, the plots shows the 
relative difference between the numerically computed and $ for a single Kerr solution with the same mass and 
angular momentum. The initial location of the horizon is indicated. These plots demonstrate that outside the horizon 
$ approaches the Kerr form as the separation between the throats (i.e., between the singularities in $sing) is reduced. 
The solutions deep inside the horizon are very different, as they must be, no matter how small the separation. As 
explained above, for these choices there will be a formal singularity at fz corresponding to a TZ^ x feature. 

So far we have presented initial data solutions in this section which can be used for close-limit analytic work, and 
in Sec. II we presented solutions that accurately describe two Kerr holes in the far separation limit. In the remainder 
of this section we will show that, with very little additional complexity, these close and far limit solutions can be 
combined. To simplify the description of how this is done, we will confine attention to the case of two identical holes, 
of parameter Mi = M2 = M/2 and angular momentum Ji = J2 = J/2, symmetrically placed at zi = —Z2 = zq, as in 
the examples of Sec. I. 

To combine solutions we need to introduce a mixing function /(zq/M) that is used to weigh the amount of the 
far-limit solution that is to be mixed with the close limit solution. For the choice made for Ji/Mi, let us take Zcrit 
to be the largest value of zq for which the initial horizon encloses all throats (the two throats representing the holes, 
and the feature at r = 0). Then the mixing function must have the property that / = 1 for zq/M > Zcrit/M. It must 
also have the close-limit property that / — > for zq/M — > 0. The mixing function might be taken, for example, to be 

. f 1 for Zq > Zcrit 

■ ~ \ 1 - (Zo - Zcritf lzl„i. for Zo < Zcrit 

For our mixed solutions, we take q and the H functions to be 

q = /[ft (Ml, Ji) + g2(Af2, J2)] + (!-/) qK{M, J) , (28) 
He = /fci(Mi, Ji) + He2{M2, J2)] 

+ f)HEK{M,J) , (29) 
Hp - f[Hpi{Mi, Ji) + Hp2iM2, J2)] 

+ {1- f)HpK{M,J) . (30) 



For clarity we have written the expressions in (g8|) - (|30|) with parameters Mi , Ji , M2 , J2 , M, J, and have not invoked 
the special values of Mi, M2, Ji,J2- The functions (gi, Hpi, Hpi), with subscript 1, are the shifted solutions of (p^), 
(pT|), and (p^, and similarly for the functions with subscript 2. The functions with K subscript are the "pure" Kerr 
functions from (H)-(^). For the mixed solution we take the function $sing to have the form 
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$si„g = 1+2 



1 f^/W 



aV2 1 /y/Ml 



- + (1 _ 1)7^2 



— a- 



;V2 



(31) 



•y/r^ — 2zir cos 9 + zf 



2 



\Jr'^ — 2z2r cos + z| 



For separations large enough so that a horizon does not enclose both holes, the solution, is identical to the solution 
of Sec. I, and therefore has the correct far limit, i.e., in the limit of large separation the solution near each throat 
approaches that of an isolated Kerr throat. For zt^jM <C 1, the mixing function / approaches 0, and the solution 
approaches the Kerr solution, and hence the mixed solution must have the correct close-limit, i.e., in the limit of small 
separation the solution outside the horizon approaches that of a single Kerr hole. 

It is interesting to note that our sequence of solutions, with mass and the angular momentum the same at large 
and at small separations, cannot produce a rapidly rotating final hole; for extreme individual holes, with M\ = Ji, 
the final hole has a = J /M"^ = 0.5. If we want a close limit final hole that is rapidly rotating we cannot use a 
constant-mass, constant-angular momentum sequence. This reminds us that simply adding rotating holes does not 
automatically produce something of physical relevance, or at least, does not produce what we may want. This is a 
point we will be coming back to in the next section. 



We have presented a general method for constructing axially symmetric superpositions of Kerr black holes. The 
crucial feature of this method is the adoption of a 3-metric form, that of with which the momentum constraint 
for the conformally related extrinsic curvature reduces to a linear equation, independent of the 3-geometry. In this 
limited sense, the approach has the simplicity of the Bowen-York |^ conformally flat approach. 

In the case of sizeable separation between the superposed holes, the method produces holes with "local Kerricity," 
and has the advantage of great simplicity and easy numerical implementation. The method has been illustrated with 
numerical results for two equivalent (equal mass, equal spin) holes, but many extensions are possible. It should be 
clear how to superpose any number of holes with "local Kerricity," or holes without local Kerricity, in particular 
locally Bowen-York holes. 

A crucial feature of our superposed solutions is that their close-limit is not a single Kerr hole outside the horizon. 
We have argued that unlike the spherically symmetric case, it should be expected on physical grounds that two locally 
Kerr holes, no matter how close they are, will not superpose into a single Kerr hole. Though the "failure" of the close 
limit is physically correct, it is inconvenient. We would like to do close limit perturbation theory, which will at the 
very least be of value for checking numerical codes for colliding rotating holes. We wish to emphasize that the correct 
close limit data for two holes, is the data produced by the evolution of a pair of holes falling into each other. A guess 
at what that initial data solution is, has little chance of being accurate. This was true, of course, also for initial data 
solutions representing nonrotating black holes, but in that case the variety of choices was small, and the radiation 
evolving from the close limit initial data seemed to be rather insensitive to the choice of details. This is very unlikely 
to be the case for rotating holes; here the richness of the geometry and extrinsic curvature gives a much greater scope 
of choices. 

For the reasons just given, our focus is on close limit calculations as a check on numerical relativity calculations, so 
definitiveness and simple implementation are of importance. We have given, in Sec. Ill, an alternative superposition 
scheme, which is very easily implemented, and which has the correct close limit, but it has a singularity at the 
midpoint between two equal holes. That singularity will be hidden deep inside the horizon of a close limit solution, 
and is of concern only for numerical relativity methods that require a solution throughout the interior of the horizon. 
The presence of the singularity does not in any way suggest that the solution outside the horizon is a good or bad 
representation of the late stage of infall of two Kerr holes, or whether the radiation produced outside the horizon is 
for any reason a good or bad test of numerical relativity. 

Having the singularity hidden inside the horizon means of course, that this type of initial data is only applicable in 
the close limit. Since it may be useful to have a sequence of solutions which includes both large and small separations, 
we have also given, in Sec. Ill, a method of mixing both types of superpositions so that for large separations one has 
locally Kerr holes and no spurious singularities, and in the close limit one approaches a single Kerr hole, and has a 
spurious singularity hidden inside the horizon. 

We intend to use the close limit, and mixed-sequence, initial data solutions as a starting point for close limit 
computations of radiation. We have already converted the superposition solutions, in the close limit, to Cauchy data 
for the Teukolsky equation |Q and will be evolving a number of examples. With the technical tools in place to do 
studies of evolution to a rapidly rotating hole, we will be able to answer a number of interesting questions, such as 
the range of validity of the slow rotation limit. 



IV. SUMMARY AND CONCLUSIONS 
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FIG. 1. (a) Numerically computed $rcg" for one single Kerr black hole with M = 0.5, a = 0.45, located at f = zq = 0.1 
and = 0. Details of the functional behavior near the singularity are displayed in the inset, (b) The relative difference, 
A = ($reg™ — ^reg)/$reg, between the numerical result and the analytic expression for the regular part 3>reg, for a single 

shifted hole. The difference is is shown in a neighborhood of the singularity. 
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1.6t 




FIG. 2. (a) $r™, computed for two black holes with masses Mi — M2 = 0.5 and angular momenta Ji = J2 = 0.225, 
separated by a coordinate distance f = 2zo = 4.0. (b) The relative difference, A = ($reg™' — 4>reg)/$reg, between the computed 
function 'E'"™ and the analytic expression for 3>reg for a single shifted Kerr hole of the same mass and angular momentum. See 
text for details. 
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FIG. 4. The conformal factor in the close limit. Part (a) shows cj>""™^ the conformal factor computed for two holes with 
equal masses Mi — M2 — 0.5 and equal angular momenta Ji = J2 ~ 0.225, placed symmetrically at Zi = 0.02 and Z2 ~ —0.02. 
Part (b) shows the fractional difference A = ($""™ — ^k)/^k between the two hole solution and the conformal factor $k for 
a Kerr hole of the same angular momentum and mass. See text for details. 
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FIG. 5. The fractional difference A = ($"™ — ^k)/^k between the numerically computed conformal factor for 

M = 1.0, J = 0.45 and ^k, the conformal factor for a Kerr hole of the same mass and angular momentum. In (a), the separation 
is given by 2zo = 0.2, and the difference is less than 8% outside the horizon (shown as a thick curve). In (b), 2zo = 0.8, and 
there is a sizeable difference between the close-limit and the Kerr conformal factors. 
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